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ABSTRACT 

We present measurements of the specific ultraviolet luminosity density from a sample of 483 galaxies at 
6 < z < 8. These galaxies were selected from new deep near-infrared Hubble Space Telescope imaging from 
the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey, Hubble UltraDeep Field 2009 and 
WFC3 Early Release Science programs. In contrast to the majority of previous analyses, which assume that the 
distribution of galaxy ultraviolet (UV) luminosities follows a Schechter distribution, and that the distribution 
continues to luminosities far below our observable limit, we investigate the contribution to reionization from 
galaxies which we can observe, free from these assumptions. Due to our larger survey volume, wider wave- 
length coverage, and updated assumptions about the clumping of gas in the intergalactic medium (IGM), we 
find that the observable population of galaxies can sustain a fully reionized IGM at z = 6, if the average ionizing 
photon escape fraction (f esc ) is ^30%. A number of previous studies have measured UV luminosity densities 
at these redshifts that vary by a factor of 5, with many concluding that galaxies could not complete reionization 
by z = 6 unless a large population of galaxies fainter than the detection limit were invoked, or extremely high 
values of f esc were present. The observed UV luminosity density from our observed galaxy samples at z = 7 
and 8 is not sufficient to maintain a fully reionized IGM unless f esc > 50%. We examine the contribution from 
galaxies in different luminosity ranges, and find that the sub-L* galaxies we detect are stronger contributors to 
the ionizing photon budget than the L > L* population, unless f esc is luminosity dependent. Combining our 
observations with constraints on the emission rate of ionizing photons from Lya forest observations at z = 6, we 
find that we can constrain f esc < 34% (2a) if the observed galaxies are the only contributors to reionization, or 
< 13% (2a) if the luminosity function extends to a limiting magnitude of Mr/y = -13. These escape fractions 
are sufficient to complete reionization by z = 6. Current constraints on the high-redshift galaxy population 
imply that the volume ionized fraction of the IGM, while consistent with unity at z < 6, appears to drop at red- 
shifts not much higher than 7, consistent with a number of complementary reionization probes. If faint galaxies 
dominated the ionizing photon budget at z = 6-7, future extremely deep observations with the James Webb 
Space Telescope will probe deep enough to see them, providing an indirect constraint on the global ionizing 
photon escape fraction. 

Subject headings: early universe — galaxies: evolution — galaxies: formation — galaxies: high-redshift — 
ultraviolet: galaxies 



1. INTRODUCTION 

The reionization of the intergalactic medium (IGM) is the 
last major phase transition in the Universe, and as such has 
been a major topic of recent study. Over the past decade, 
a number of lines of evidence have yielded improved con- 
straints on the duration of reionization, the completion red- 
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shift, and the primary source of ionizing photons. 

Observations of the Cosmic Microwave Background 
(CMB) measure the Thomson scattering optical depth due 
to electrons in the path. The updated 7 year Wilkinson Mi- 
crowave Anisotropy Probe results measure r = 0.088 ± 0.014, 
which corresponds to a redshift of insta ntaneous reionization 
of z = 10.6 ± 1.2 (iKomatsu et al.ll2011l) . However, reioniza- 
tion is likely to be a much more extended event, especially 
if faint galaxies play a dominant role. Measurements of the 
kinetic Sunyaev-Zeld ovich (kSZ) effect with the South Pole 
Telescope (SPT) from lZahn et al.l (1201 ll) have recently placed 
a limit on the duration of reionization, of Az < 7.9. When 
combined with the Wilkinson Microwave Anisotropy Probe 
(WMAP) constrains on instantaneous reionization, they con- 
clude at 95% confidence that reionization was complete at 
z > 5.8. Additionally, observations of the Lya forest in high- 
redshift quasars also imply an end to reionization at z ~ 6, 
as the measurements of the near-zones around z = 6 quasars 
yield results co nsistent with a negligible neutral fraction (e.g., 
iFan et al.1l20 06), though more z > 6 quas ar sight-lines ma y be 
nee ded to conclusively verify this result (M esingerl2010l) (see 
also iBecker et al.ll2007l) . 

Nonetheless, some sources must be responsible for reion- 
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izing the IGM between 6 < z < 14. Quasars were a nat- 
ural choice, as they are extremely bright, and the bulk of 
their ionizing photons can escape. However, the quasar 
luminosity function peaks at z ~ 2, and falls of f rapidly 
toward higher redshift (e.g., Hopkins et al.l 12007b . Addi- 
tionally, observations of the X-ray background rule out a 
dominan t contribution to the r eionizing photon budget from 
quasars (iDijkstra et al.l 1200 4a). Po pulation III stars, due to 
their predicte d high masses (e.g., iBromm & Larsonl 120041 
lGlover || 2005l). were very efficient emitters of ionizing pho- 
tons (ITumlinson & Shullll2000t IBromm et al.ll200lt ISchaererl 
120021) . Their overall contribut ion to reionization may , how- 
ever, have been limited (e.g.. iGreif & Bromml 120061) . This 
conclusion has been strengthened by recent simulation results 
which have corrected th e mass scale of the first stars down 
to les s extreme values dClark et al.l 120111; IGreif et al.l 1201 1L 
l20T2l) . 

The most likely source for the bulk of ionizing photons 
is thus from star-forming galaxies themselves, where ioniz- 
ing photons are created from the massive stars present dur- 
ing ongoing star formation. This has been well-studied ob- 
servationally, but it was only in the last decade when large 
samples of galaxies close to the reionization epoch were com- 
piled. Following wide and deep surveys with the Advanced 
Camera for Surveys (ACS) onboard the Hubble Space Tele- 
scope (HS T), large samples of z ~ 6 galaxies wer e com- 
piled (e.g..lGiavajisco etai]|2004tlBunker et al.ll2004l) . While 
iGiavalisco et al.l (120041) found that the specific UV luminos- 
ity density (and thus the star-f ormation rate densit y) was 
roughly constant from z = 4 to 6, iBunker et al.l (120041) found, 
via deeper observations, that this quantity appeared to de- 
cline towards higher redshift. Using the commonly ac- 
cepted (at the time) large value for the clumping factor of 
the IGM, they concluded that galaxies were unable to ac- 
count for the necessary ionizing photons to reionize the uni- 
verse. lYan & Windhorst! (120041) pointed out that if the faint- 
end slope of the galaxy luminosity function was steep, it could 
be that L < 0. IL* galaxies dominated the ionizing photon bud- 
get, and thus allowed galaxies to complete reionization by z = 
6. Regardless of the interpretation, these galaxy samples con- 
tained only ~ 50-100 galaxies, each of which was only de- 
tected in one band, given the lack of deep HST near-infrared 
observations at the time. 

The advent of the Wide Field Camera 3 (WFC3) on HST 
has opened the door to the z > 6 universe, allowing the first 
investigation into whether reionization completed at z = 7 or 
earlier. This is combined with recent results which have vastly 
changed the expected IGM clum ping factor to much lower 
values. [Finkelstein et al.l (120101) investigated the combined 
UV luminosities of their observed z = 7 galaxies, and found 
that while they could come close to reionizing the IGM, either 
fainter galaxies or high (> 50%) escape fraction s were neces- 
sary t o complete reionization by z = 7 (see also lBunker et al.1 
l2010h 

A number of studies have used measured luminosity func- 
tions of z ~ 7 and 8 galaxies to infer their contribution to 
reionization, by integrating them down to an assumed mag- 
nitude limit. The majority of these studies found that, once 
galaxies fainter than the detection threshold are accounted 
for, galaxies could reionize the universe at such early times 



(e.g. lOesch et al.l20105lMcLure et al.l2010tlLorenzoni etaTl 



1201 U IBouwens et all 1201 lbl) (see also iGrazian et all 1201 ll) . 
This is however reliant on a number of assumptions. First, 
the luminosity function results are susceptible to the assump- 



tion of a Schechter function parameterization. At z ~ 5- 
6, where the samples are larger, there is no strong evidence 
for deviation from this function (e.g.. iBouwens et all 120071: 
iMcLur e et al. 20091). However, as we push closer to the Big 
Bang, galaxies are changing rapidly, thus at some point we 
may encounter an epoch where the Schechter function is no 
longer an accurate representation. This could be due to a va- 
riety of effects, with one being a lack of active galactic nuclei 
(AGN) feedback if AGNs are not yet present in the centers of 
all galaxies at very high redshifts, which is plausible depend- 
ing on the speed with which supermassive black holes form 
in the early universe. 

Additionally, while recent evidenc e has indicated a very 
steep faint end sl o pe at z > 7 (e.g.. IBouwens et al J 1201 let 
lOesch et all 120121: Bradley et ail 120121) . the uncertainty on 
these measurement s are large (e.g., <r(a) sa 0.2 at z = 8 
IBradlev et al.ll2012h . Finally, when integrating the luminosity 
function, one needs to choose a limiting magnitude, as sup- 
pression from the UV background will result in the gas in 
galaxies becoming heated below some circular velocity (i.e. 
mass) limit. Galaxies above this limit are dense enough where 
collisional excitation of H I by electrons can overcome pho- 
toionization heating, and have their gas cool and form stars. 
However, we have no observational evidence for this value, 
and theoretical res ults from the literature yield values of -15 
<M Um < -10 (e.g..lFlnlator et al.l201 UlMufloz & Loebl201 U 
iKulkarni & Choudhury|201 UlChoudhurv et al.ll2008l ). With a 
steep faint-end slope, a difference of this level can result in a 
difference in the integrated luminosity density by more than a 
factor of two. 

Here we measure the specific luminosity density of the ob- 
servable galaxy population using the largest sample of 6 < 
z < 8 galaxies yet compiled, which allows us to study the con- 
tribution of galaxies to reionization without the uncertainties 
inherent in invoking a parametrized luminosity function. This 
allows us to assess whether a significant contribution from 
galaxies below the observational limits is necessary to com- 
plete reionization at a given redshift. We emphasize our re- 
sults at z = 6, where we have a large sample, and each galaxy 
is detected in four individual imaging bands, yielding much 
more robust results over the previous studies which had ac- 
cess only to optical ACS data. 

In §2 we describe the datasets used in this study, as well 
as our photometry and sample selection methods. We dis- 
cuss how we measured the rest-frame UV specific luminos- 
ity density, and corrected it for incompleteness down to our 
magnitude limit in §3, while in §4, we discuss the implica- 
tions of our results on the reionization of the IGM. Through- 
out this paper we assume a concordance cosmology, with H 
= 70 km s" 1 Mpc" 1 , Q m = 0.3 and fl x = 0.7. All magni- 
tudes are reporte d in the AB system , where m\s = 31.4 — 
2.51og(/„/lnJy) dOke & Gunnlll983l) . In the remainder of 
the paper we make use of luminosity functions from the liter- 
ature over the redshif t range of our s t udy. W e use the lumi- 
n osity functions f r om IBouwens et all (120071) for z = 4, 5 and 
6; IBouwens etail (1201 Id) for z = 7; IBradlev et al.l (120121) for 
z = 8 and IBouwens et al.l (1201 lal) for z = 10. These represent 
the most up-to-date luminosity functions for these redshifts at 
the time of this writing. The characteristic magnitude values 
from these studies are: M* = -20.98 (z = 4), -20.64 (z = 5), 
-20.24 (z = 6), -20. 14 (z = 7), -20.26 (z = 8) and -18.3 (z = 
10). 
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2. OBSERVATIONS 

In this Paper, we use a sample of high-redshift galaxies se- 
lected from HST data in the GOODS-S region. We review our 
methods below, and not e that they are similar to those used in 
iFinkelstein et al.l d201 lb . thus the reader is referred there for 
further details. 

2.1. Data 

We use WFC3 imaging data from three separate programs. 
The first is the HUDF09 program (PID 11563; PI Ming- 
worth), which obtained deep imaging over the Hubble Ultra- 
Deep Field (HUDF) as well as its two flanking fields (which 
we will refer to as HUDF09-01 or PARI and HUDF09-02 
or PAR2). The second dataset comes from the WFC3 Early 
Release Science Prog ram (ERS; PID 11359; PI O'Connell; 
IWindhorst et al.ll201 ll) . which imaged the northern ^25% of 
the GOODS-S field. The third dataset comes from the Cos- 
mic Assembly Near-Infrared Deep Extra galactic Legacy Sur- 
vey (CANDELS; PI F aber & Ferguson HGrogin et aljfeoi U 
iKoekemoer et al.|[201 lb . which obtained deep imaging of the 
central ~50% of GOODS-S (hereafter referred to as CAN- 
DELS DEEP), and less deep imaging of the southern ~25% 
of GOODS-S (hereafter referred to as CANDELS WIDE; PID 
12060, 12061 & 12062). We use the completed CANDELS 
dataset in GOODS-S, which includes 10 epochs in the DEEP 
field, as well as both epochs of WIDE imaging. All three of 
these surveys obtained imaging in the F125W and F160W fil- 
ters (hereafter J 125 and Hi 6 o). The CANDELS and HUDF09 
surveys also obtained imaging in the F105W filter, while the 
ERS used the narrower F098M filter (hereafter Y 105 and Y098, 
respectively). These fields and their associated depths are 
summarized in Table 1 . 

2.2. Photometry and Sample Selection 

We performed photometry on all three datasets us ing the 
Source Extractor software (iBertin & Arnoutsl 11996b . We 
chose a single set of parameters that we applied to all im- 
ages, which we tested to maximize detection of real sources 
while minimizing contamination (the key parameters are 
DETECT_MINAREA=7 and DETECT_THRESH=0.6). We 
did this in two image mode, using a weighted sum of the 
J125 and H160 images as the detection image, and the three 
WFC3 bands, as well as the archival GOODS-S ACS F435W, 
F606W, F775W and F850LP (hereafter B435, V 60 6, i?75 and 
zs5o, respectively) bands as the measurement images. We also 
measured photometry on very deep F8 14W ACS imaging, ob- 
tained in parallel during the CANDELS observations. Re- 
cently obtained ACS imaging suffers from poor charge trans- 
fer efficiency; this can be corrected for detectable sources 
dAnderson & Bedinl 120101) . but at this point it is unclear 
whether robust upper limits can be obtained from these data. 
We thus use the F814W imaging (which has had this correc- 
tion applied) only as a veto band when inspecting z = 8 galaxy 
candidates, as at these redshifts there should be no detectable 
flux in this band. All images were matched to the point-spread 
function (PSF) of the H160 image prior to the photometry. 

We selected our sample of high-redshift ga laxies using the 
photo metric redshift fitting software EAZY (iBrammer et al.l 
2008). Rather than using the best-fit photometric redshift, we 
used the full probability distribution function (PDF) to select 
our galaxies. We do this by compiling samples in three red- 
shift bins, z = 6, 7 and 8. To make it into a given bin, an object 
must have a > 3.5er detection in both the J125 and H160 bands, 
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FIG. 1 . — The number of objects selected in each of our six fields for z = 6, 
7 and 8, from left-to-right. At z = 6 - 7, the CANDELS Deep and ERS fields 
dominate the sample. However, at z = 8, due to the decrease in the char- 
acteristic magnitude M* and the increased luminosity distance, the deeper 
HUDF09 fields dominate the sample. The CANDELS Wide field contributes 
a non-negligible number of galaxies at z = 6, but is not deep enough to con- 
tribute significantly at higher redshifts. 

it must have > 70% of its integrated PDF in the primary red- 
shift solution, and it must have >25% of its integrated PDF 
within z S ampie± 0.8 (i.e., for z = 6, >25% of the integrated 
PDF must be at 5.2 < z < 6.8). Finally, objects which satisfy 
this criterion for multiple samples are placed in the sample 
which contains a larger fraction of the PDF (for this purpose, 
we also integrated the PDF for hypothetical z = 5 and z = 9 
samples, to prevent those objects from being placed in our z = 
6 and z = 8 samp les, respectively). These methods are simi- 
lar to that used in Fink elstein et al.l (1201 lb . and their Figure 1 
illustrates this selection technique. 

This selection was performed separately for each imaging 
dataset. All selected candidates were visually inspected, to 
screen against false detections such as diffraction spikes or 
oversplit regions on the edges of bright galaxies. Additionally, 
we screened against stellar contaminants, first by investigat- 
ing whether a given source is resolved, which we did by com- 
paring the measured full-width at half-maximum (FWHM) of 
the sources to that of stars in the image (where this measure- 
ment was performed by Source Extractor). For sources that 
had FWHMs close to that of stars in the image, we compared 
their colors to measured colors of late-types stars and brown 
dwarfs, integrating observed NASA Infrared Telescope Facil- 
ity (IRTF) SpeX prism spectra through the ACS and WFC3 
bandpasses 1 . We found 15 likely M-dwarfs in our z = 6 sam- 
ple, and 4 likely brown dwarfs in our z = 7 sample. After 
clearing these contaminants from the sample, we had a total 
of 302, 136 and 45 galaxy candidates at z = 6, 7 and 8, re- 
spectively. These samples are summarized in Figure 1, which 
highlights the importance of a wide area at z — 6, and very 
deep imaging at z = 8. 

3. ANALYSIS 

3.1. Measuring the Rest-Frame UV Light 

In order to study the contribution of the galaxies in our sam- 
ple to the reionization of the universe, we need to measure 
their rest-frame UV absolute magn itudes. We have chosen 
to do this in the same manner as in IFinkelstein et al.l (1201 lb . 
which we briefly summarize here. We first compute a grid 
of synthetic stellar population models over a range of red- 
shifts, population ages, stellar metallicities, star-formation 
histories and dust a ttenua tions, using the updated models of 
iBruzual & Charlotl (120031 hereafter BC03). We then found 
the best fit model for a given galaxy in our sample via \ 2 

1 This research has benefitted from the SpeX Prism 
Spectral Libraries, maintained by Adam Burgasser at 
http://pono.ucsd.edu/~adam/browndwarfs/spexprism 
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TABLE 1 
Observations Summary 



Field 


Area 


B435 


V606 


1775 


1814 


z 850 


^098 


Y U >5 


^125 


#160 




(arcmirr ) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 



HUDF09 MAIN 


5.0 


29.5 


29.9 


29.7 




29.0 




29.2 


29.5 


29.5 


HUDF09 PARI 


4.4 




28.9 


28.7 




28.6 




28.0 


28.9 


28.7 


HUDF09 PAR2 


4.6 




29.0 


28.7 




28.4 




28.7 


29.1 


28.9 


GOODS-S DEEP 


62.0 


28.1 


28.3 


27.7 


28.2 


27.5 




28.2 


28.1 


27.9 


GOODS-S WIDE 


34.1 


28.1 


28.3 


27.7 


27.8 


27.5 




27.4 


27.6 


27.3 


GOODS-S ERS 


40.5 


28.1 


28.3 


27.7 


27.7 


27.5 


27.6 




28.0 


27.7 



NOTE. — The area is measured from the image area used to detect objects, which excludes noisy 
regions on the edges, regions without data from all available filters, as well as regions overlapping 
deeper data (i.e., the CANDELS Deep field also covers the HUDF, thus we have masked out this re- 
gion in the CANDELS imaging for this analysis). The remaining columns are 5a limiting magnitudes 
measured in a 0.4"-diameter apertures on non-PSF matched images. The F814W data from CAN- 
DELS was only used in the visual inspection of z = 8 candidate galaxies, primarily in the GOODS-S 
DEEP field as it is ~0.5 mag deeper than the existing GOODS ACS data. 
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z z z 

FIG. 2. — A summary of our completeness simulations. Each panel represents one of our redshift bins of z = 6, 7 and 8 from left-to-right. Each panel contains 
the probability of recovering a given galaxy (as a function of redshift and M1500) for two magnitude values; Misoo = —22 (solid line), and M1500 = —20 (dashed 
line). At z = 6, 7 and 8 M1500 = —22 corresponds to apparent UV magnitudes of 24.7, 24.9 and 25.1, and M1500 = —20 corresponds to 26.7, 26.9 and 27.1, 
respectively. While all fields can efficiently detect bright galaxies at any redshift, only the two HUDF09 fields can efficiently detect fainter galaxies at higher 
redshift, which is consistent with our recovered distribution shown in Figure 1 . Although the HUDF09 fields are more efficient than the CANDELS Deep and 
ERS fields for fainter galaxies at z = 6, the much larger area probed by the latter fields results in their dominating the sample at that redshift. We note that although 
they are of comparable area, the volume probed by the WIDE field is much less than the ERS. At fainter magnitudes, this is due to the increased depth in the 
ERS. However, even at bright magnitudes, the ERS is more complete at z = 7 and 8, as the F098M filter used in the ERS does not overlap with F125W (while 
F105W, used in the remaining fields, does), allowing more precise photometric redshifts. 



minimization, and then integrated this best-fit model spectrum 
through a 100 A-wide square bandpass centered at rest-frame 
1500 A to compute the flux at that wavelength. This was then 
converted to the absolute magnitude at 1500 A (M1500) using 
the distance modulus to the best-fit photometric redshift. The 
uncertainty on M1500 was derived via a series of 100 Monte 
Carlo simulations, where in each simulation each galaxy's 
fluxes were altered by a Gaussian random distribution with a 
mean of zero and a standard deviation equal to the photomet- 
ric error for the band in question. In each simulation, a new 
redshift was chosen based on the photometric redshift PDF for 
a given object; in this way the uncertainties on all derived pa- 
rameters include the inherent uncertainty in the redshift. The 
68% confidence range of M1500 was then derived as the cen- 
tral 68% spread of M 15 qo values from the simulations for each 



object. 



3.2. Incompleteness Correction 



Prior to computing the specific luminosity density, we first 
need to correct for incompleteness down to our limiting mag- 
nitude. We have done this by inserting mock galaxies into 
our images, running separate simulations for each redshift bin 
and each dataset (for a total of 18 sets of simulations). In 
each simulation, we chose a uniform random distribution of 
redshifts ranging from z sa mple ± 1, as well as a range of H160 
magnitudes modeled after the number counts in a given field, 
but continuing the trend to fainter magnitudes. We then chose 
a stellar metallicity, population age, and color excess E(B-V) 
from distributions which were tuned such that the recovered 
galaxies in a given simulation had a similar rest-frame UV 
spectral slope distribution as those of the recovered galaxy 
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sample (ensuring that we're not correcting for, e.g., a pop- 
ulation of red galaxies that do not exist in reality). A stellar 
population model was then created with these parameters, and 
integrated through the bandpasses of interest to derive galaxy 
magnitudes in each band. 

The image of each m ock galaxy was c reated with the GAL- 
FIT software package dPeng et al.ll2002l) . The position angle 
was drawn from a uniform random distribution, while the Ser- 
sic index n was drawn from a log-normal distribution, such 
that the mock sample was predominantly disk galaxies, with 
fewer bulge-dominated galaxies. Both the half-light radius 
and the axial ratio were tuned such that the recovered distri- 
butions matched the distributions from our observed galaxy 
sample (to prevent against, i.e., correcting for large galaxies 
which do not exist). 

Samples of 100-200 mock galaxies were placed in sub- 
regions of the images in each field. The images were then 
treated in the same manner as in §2.2, in that we derived 
photometry with Source Extractor, and photometric redshifts 
were measured with EAZY. Samples were selected at each 
redshift using the same criteria as that used on our galaxy sam- 
ples. This process was repeated until ~ 10 5 galaxies were in- 
put and measured from the images (i.e., 500-1000 iterations). 
Typical numbers of recovered galaxies were ~ 2 - 3 x 10 4 
(depending on the field and redshift), where a galaxy is de- 
noted as recovered if it is detected in the photometry and se- 
lected as a candidate in the redshift range of interest. 

For each field and at each redshift, we then computed the 
probability that a given galaxy was recovered in our sample as 
a function of the input redshift and Mi 500 (where here Mi 500 
was measured from the input spectrum to the simulations), 
defined as V(M,z). These values are shown at a few represen- 
tative magnitudes in Figure 2. 

3.3. Evolution of the Rest- Frame UV Specific Luminosity 
Density 

Using the value of M1500 for each galaxy, as well as V(M,z) 
for each field and redshift, we computed the rest-frame UV 
specific luminosity density (puv), m units of erg s -1 Hz" 1 
Mpc" 3 (specific refers to per unit frequency, while density 
refers to per unit volume). We converted the absolute mag- 
nitude M1500 of each galaxy to the 1500 A specific luminos- 
ity (i.e., erg s -1 Hz" 1 ) using the best-fit photometric redshift 
(while noting that the uncertainty in M1500 includes the uncer- 
tainty in the redshift; see §3.1). 

We computed the effective volume, which is the volume 
over which we are selecting galaxies given our imaging depths 
and selection criteria, defined as 

fdV 

V ef f(M l5 oo)fieid= Af ie i d / —V(M l5Q0 ,z)dz (1) 

where A is the area probed by a given field, which is given 
in Table 1, and ^ is the co-moving volume per square ar- 
cminute per redshift slice dz at the redshift z- The value of 
puv is thus the sum of the specific luminosities divided by 
the effective volume for a given field and redshift, in each 
magnitude bin. We note that in the following, when we refer 
to measurements of our observed galaxy sample, or the ob- 
servable galaxy sample at a given redshift, we refer to these 
incompleteness-corrected measurements. 

We computed a total value of puv for each redshift, by first 
taking the sum of all specific luminosities in a given redshift 



sample in a given magnitude bin, and then dividing it by the 
sum of the effective volumes for each field in that magnitude 
bin. We corrected down to a limiting M1500 = -18, which cor- 
responds to 0.13, 0.14 and 0.1 9 x L* atz= 6, 7 an d 8, respec- 
tively (as well as 0.06 x L* 3 : ISteidel eTaT1(fT99l '). The total 
value of puv was then the sum of these values over all magni- 
tude bins. These values are shown as the large cyan circles in 
Figure 3, and are tabulated in Table 2. 

The uncertainties quoted on pjjy are comprised of three 
components: Poisson noise, photometric scatter, and uncer- 
tainty in the effective volumes. These were included simulta- 
neously in a series of 10 4 bootstrap Monte Carlo simulations. 
In each simulation, we accounted for Poisson noise by select- 
ing A^' = N + A x \/N galaxies in each redshift bin through 
random sampling with replacement, where N is the number 
of galaxies in the bin, and A is a random number drawn from 
a normal distribution with a mean of zero and a standard de- 
viation of unity. We then account for the photometric noise 
by taking each of the galaxies in the modified sample, and 
obtaining a new estimate, M( 500 , for M1500, where M' 15O0 is 
the original value Muoo modified by a Gaussian with a equal 
to the photometric error (performed in flux-space). Finally, 
we accounted for the uncertainty in the effective volume sim- 
ulations by computing the uncertainty in V(M\^^z) via a 
separate set of simulations, randomly selecting N recovered 
simulated galaxies, where N is 50% of the total number of 
recovered galaxies. This was done 10 3 times, generating 10 3 
■P(Mi5oo,z) arrays. The uncertainty in 'P(Mi5oo,z) was derived 
as the standard deviation of the 10 3 probability values at each 
magnitude and redshift. These uncertainties were then folded 
in to our uncertainty on puv during the Monte Carlo simu- 
lations by modifying each value of the ^(Misocz) array by a 
Gaussian with a equal to the uncertainty in T^M^oo^)- Thus, 
at each redshift we obtained 10 4 values for puv including the 
effects of sample size, photometric scatter, and uncertainty in 
the effective volumes. The uncertainty we quote is the central 
68% range of these values. 

We list the total effective volumes probed by our study at 
each redshift (at two representative magnitudes) in Table 2. 
We note that these volumes are subject to our assumptions of 
the galaxy properties in our incompleteness simulations. As 
discussed in the previous subsection, the key parameters in 
recovering a galaxy are the color and size. While we have 
chosen what we believe are appropriate distributions of both 
quantities (see §3.2 for discussion), it is however interesting 
to examine how the effective volume would change if the as- 
sumptions were different. We have thus re-run our simula- 
tions at z = 6 in each field, using a modified color and size 
distribution (in comparison to our default distributions, we 
broadened the distribution of dust attenuation, and changed 
the size distribution to peak at smaller sizes, but with a higher 
large-radius tail). We found that the derived effective volumes 
from this modified simulation were ~ 10% smaller than our 
nominal simulations in §3.2. Thus, the derived UV luminosity 
densities would be ~ 10% higher if these volumes were more 
physically correct. However, as discussed above, we believe 
our assumed physical property distributions are physically ro- 
bust, so our effective volumes should be representative of the 
true volumes probed by our survey. As this test shows, if for 
whatever reason this is not the case, it creates a systematic un- 
certainty on puv on the order of 10% for modest changes in 
the input distributions. 
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TABLE 2 

Rest-Frame UV Specific Luminosity Density 



z 


Ngo/ cade s 


puv 


puv 


puv 


puv 


V ef f (M 1500 = -22) 


Ve//(M 15 oo=-20) 






M1500 <-18 


L > 0.2L* 


L > 0.5L* 


L > L* 


(MpcT 3 ) 


(MpcT 3 ) 


6 


302 


0.79 ± 0.06 


0.67 ± 0.05 


0.44 ± 0.04 


0.25 ± 0.03 


3.07 ± 0.06 x 10 5 


2.28 ± 0.03 x 10 5 


7 


136 


0.44 ± 0.04 


0.43 ± 0.04 


0.27 ± 0.03 


0.17 ±0.02 


3.04 ± 0.04 x 10 5 


2.23 ± 0.02 x 10 5 


8 


45 


0.27 ± 0.06 


0.27 ± 0.06 


0.14 ±0.03 


0.07 ± 0.02 


2.96 ± 0.05 x 10 5 


1.94 ±0.02 x 10 5 



NOTE. — All specific luminosity density values are in the units of 10 erg s Hz Mpc . The effective 
volumes are quoted in comoving units, and were derived from the completeness simulations discussed in 
§3.2. The uncertainties on the volumes come from bootstrap simulations, succesively selecting 50% of the 
recovered simulated galaxies to assess the variation in V(M,z)- 



3.4. Critical Specific Luminosity Density to Complete 
Reionization 

As shown in Figure 3, our results show a decrease with 
redshift in the specific luminosity density for the observable 
galaxy population down to M1500 < — 18. In order to under- 
stand the implications of this trend on the reionization of the 
IGM, we need to connect our observable, pjjy, to the neces- 
sary number of ionizing photons to sustain a fully reionized 
IGM. We do this in the sa me manner as in | Finkelstein et all 
( 120101) . where we followed iMadau et al.l (11999b to derive the 
required 1500 A luminosity density to sustain reionization: 



x Xhii — erg s" 1 Hz" 1 Mpc" 3 , 

Jesc 



(2) 



where is the cosmic baryon density and Ii-jq is the Hub- 
ble parameter in units of h = 0.7, and the constant 0.0461 is 
our assumed value of (x /i^ n ) from Wilkinson Mic rowave 
Anisotropy Probe Year 7 data dKomatsu et al.l l20llT) . The 
clumping factor of ionized hydrogen in the IGM is denoted 
by C, and ties in the average rate of recombinations in the 
IGM. The average escape fraction of ionizing photons from 
galaxies is denoted by f esc . The variable £53 is the number of 
Lyman continuum photons per unit of forming stellar mass 
in units of 10 53 photons s" 1 (M yr" 1 )" 1 . To determine this 
value, we used the updated BC03 models with a Salpeter IMF 
(with stellar masses from 0.1 - 100 M Q ) and a constant star- 
formation history. We examined values of the stellar metal- 
licity from 0.02 -1.0 Z Q finding that £53 ranges from 0.9 to 
1 A 2 . The volume ionized fraction is denoted by xhh- We start 
by assuming this is equal to unity, and investigate changes in 
§4.6. The gray curves in Figure 3 denote the critical specific 
UV luminosity density from this equation required to sustain 
a fully reionized IGM at a given redshift, where the width of 
the curve denotes the range of £53 values for differing metal- 
licities. Should the IMF be more top heavy than Salpeter, or 
the metallicities in these galaxies be <C 0.01 Z Q , it would in- 
crease the number ionizing photons per observed UV photon, 
maki ng it easier f or a given population to reionize the IGM 
(e.g.. IChary||2008l) . We note that currently there is no strong 
evidence to support such very l ow metallicities of these galax- 
ies (e.g.. iBouwens et alJl2010t iFinkelstein et alj[2010l 1201 ll 
iWilkins et al.11201 lbHFinlator et al.11201 lHWise et al.ll2012lK 

2 These values assumed an age of 50 Myr. For a constant star-formation 
history, £53 is invariant for ages larger than 10 Myr. 



We reiterate that, under the assumption of xhii = 1, these 
curves show the luminosity density required to maintain a 
fully ionized IGM. Thus, if at a given redshift our galaxies 
do not reach this threshold, it may not mean that sources of 
ionizing photons are missing, it may just indicate the epoch 
when the neutral fraction is no longer negligible. As the ion- 
ized fraction falls, one can calculate the luminosity density 
necessary to maintain a given volume ionized fraction by low- 
ering the value of xhii- Prev ious observation s constrain xhii 
to equal unity at z < 6 (e.g.. iFan et al.lr2006h . and it is likely 
negligible by z > 15 (c.f. §4.6). As we are just pushing back 
into the reionization epoch, we assume this parameter equals 



duced (possibly to ~0.5 - 0.7, e.£ 


;..lBolton & Haehneltl|2007 
l2011tlMortlocketal.l2011 


Pritchard et al.l2010tlBolton et al. 


Pentericci et al.11201 ll). The primary effect of this is that a 



lower ionized volume fraction makes it easier for galaxies to 
keep the ionized portion of the IGM in an ionized state. The 
reader can take this into account when investigating the vari- 
ous gray curves on Figure 3. 

4. DISCUSSION 

4.1. Assumptions 

The critical curves for sustaining a fully reionized IGM 
are shown as the gray bands in Figure 3, for different as- 
sumptions of the ratio between the clumping factor and 
escape fraction of ionizing photons, Clf esc . The clump- 
ing factor of ionized hydrogen in the IGM is derived the- 
oretically, and up until r ecently had been though t to be 
very large, C ~ 30 (e.g., iGnedin & Ostrikeri 119971) . This 
high value necessitates more ionizing photons to reion- 
ize a given volume, leading a number of studies to con- 
clude that galaxie s could not reionize the universe (e.g., 
iBunker et al.ll2004T) . However, recent higher-resolution sim- 
ulations yield much lower values, of C ~ 1-5, primar- 
ily due to a better understanding of the interface between 
the IGM and the circum-galactic medium, but also because 
of an impro ved understanding of the feedback from reion- 
ization (e.g iMiralda-Escude et al.l 120001: iPawlik et al.ll2009t 
Shull et al.11201 11) (see also iWise & Abell 120051: IllievetalJ 



2006tlRaicevic & Theunsl20i ll) and improvements in our un- 



derstanding of self-shielding (e.g.. lMcOuinn et afll201 ll) . Re- 
cent simulations by Finlator et al. (2012, in prep) show that 
when you only consider gas that is sufficiently under-dense 
that its ionization state is determined by the global UV back- 
ground rather than by the next nearest galaxy (the impact 
of circum-galactic gas being more appropriately attributed to 
fesc), C ~ 3 at z = 6, dropping to ~ 2 by z = 8. 
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The escape fraction of ionizing photons from galaxies is 
also poorly understood, and is one of the primary goals of a 
number of observational studies. This is impossible to directly 
observe at the redshifts of our study, as any escaping ionizing 
photons will be absorbed by lingering neutral gas in the IGM. 
Direct detection is possib le at z < 4, but the results have been 
mixed. ISiana et ail (120101) searched for ionizing photons with 
UV imaging of z ~ 1 star-forming galaxies, finding only non- 
detections, resulting in an upper limit on the relative escape 
fraction of ionizing photon s to UV ph otons of f esc .rei < 2% 
(see also ISiana et alj|2007l) . At z ~ 3, iShaplev et all ( 120061) 
found that 2/14 LBGs have relative escape fractions of ~ 50- 
100%, with the remainder having f esc ,rei < 25%. A number of 
other stu dies have detected escaping ionizing photons at z ~ 
3^1 (egJSteidel et alJl200U llwata et alj|2009l:lvanzella et alj 
120 lOt iNestor et al.11201 ll though we note some of these stud- 
ies observed in the same field (SSA22), thus they may have 
sources in common), though space-based imaging may be 
nece ssary to exclude cont amination from intervening sources 
(e.g. JVanzella et alj|2012h . 

Thus, it appears is if the average escape fraction of ionizing 
photons may be increasing with increasing redshift (or, alter- 
natively, the fraction of galaxies exhibiting f esc significantly 
larger than zero rises with redshift). Future results from UV 
observations from CANDELS, as well as in the HUDF (PI 
Teplitz) will shed more light on this issue, but there are a few 
observati onal clues to support this interpre tation. A recent 
study by iKuhlen & Faucher-Giguerel (120121) have combined 
all observational evidence (luminosity functions, WMAP, kSZ 
and Lya forest measurements) to create a concordance model 
of reionization, which implies an increasing escape fraction 
at higher redshift with f esc higher by perhaps as much as 10 x 
at z > 4. Additionally, the fraction of galaxies exhibiting Lya 
emission (with a rest-frame Lya equivalent w idth > 20 A) 
rises with redshift, from ^25 % at z ~ 3 (e.g., 
120031) to > 60% at z ~ 6 (e.g. JStark et aT1l20ia 



easily (e.g.. lKornei et al.l2012l) . In addition, the stellar masses 
of galaxies have been observed to decrease with increasing 
redshift, reducing the gra vitational potential, also allowing 
more efficient winds (e.g..lStark et a l. 2009; Finkelstein et all 
2010t IGonzalez e t al. 201^t lLabbe et alj|2010t iMcLure et all 



Shaplev et alj 
2011b. What- 



ever physical mechanism allows a larger fraction of Lya pho- 
tons to escape will likely also allow a larger fraction of ion- 
izing photons to escape (to a point; if all ionizing photons 
escape, no Lya will be produced). One possible scenario is 
that the covering fraction of neutral gas decreases with in- 
creasing redshift. Observa tional evidence for this has been 
found, as iJones et ail d201 lb measured that the low ionization 
absorption lines in a stacked rest-frame UV spectrum of z = 4 
galaxies were weaker than at z = 3, implying a lower covering 
fraction of neutral gas at higher redshift. Additionally, they 
stacked galaxies with and without Lya emission, and found 
weaker absorption features in the stack of galaxies with Lya 
emission. This supports a scenario where galaxies at higher 
redshift have a lower neutral gas covering fraction, allowing 
both more ionizing and Lya photons to escape directly to the 
observer, and that the observed increase in the incidence at 
Lya emission at higher redshift may also trace an increase in 
the ionizing escape fraction. 

The reduced covering fraction may be due to star-formation 
activity. The specific star-formation rate (sSFR; star- 
formation rat e per unit mass) is remarkably consta nt from 
z ~ 4-7 (e.g JStark et all 120091: IGonzalez et alJl20loh . How- 
ever, the typical galaxy physical size at a given luminosity 
decreases towards higher redshift (e. g., iFerguson et al. II20041 ; 
lOesch et al.ll2010atlRvan et al.ll2012b . thus the constant sSFR 
is occurring over a smaller surface area, which can result in 
the star-formation episode driving gas from the galaxies more 



201 lb . plausibly decreasing the gas covering fraction, allow- 
ing a higher escape fraction. We note however that these 
measurements solely constrain the stellar masses, not the to- 
tal baryonic mass, which may not evolve as strongly. Fi- 
nally, a higher escape fraction is also consistent with the 
observed galaxy luminosity density and the photoioniza- 
tion rate inferred from observations of the Lya forest in 
high redshift quasars, wh ich imply f esc ~ 20-30% (see §4.5 
iBolton & HaehneHll2007h . 

A number of theoretical studies have examined the escape 
of ionizing photons, both as a func t ion of redshift and halo 
mass. For example, lYanma et al.1 (1201 11) investigated this 
quantity from 3 < z < 6. At z = 6, they found that the es- 
cape fraction was > 10% for all halos more massive than 10 1 1 
M Q (which likely encompasses the majority of the galaxies 
in our sample), with f esc risin g towards lower halo masses, 
for an overall average of 40%. R azoumov & Sommer-Larsenl 
( 12010b investigate this at higher redshifts, and found that typ- 
ical escape fractions in halos with masses from 10 10 - 10 11 
M at z = 6-8 are ~ 50-80%. T he radiative transfer simu- 
lations of Haa rdt & Mada ul ( 120121) also find a relatively high 
escape fraction, and one which rises with redshift. They find 
that the luminosity weighted escape fracti on rises from 10 - 
30% fr om z = 6 — >• 8. Another recent study. IConrov & Kratterl 
(12012b investigated the impact of runaway massive stars on 
the overall escape fraction. They found that by accounting 
for the motions of these stars the escape fractions in lower 
mass galaxies could rise, as the smaller sizes of lower mass 
galaxies would allow massive stars to emerge from the gas 
disk prior to exploding. They found that galaxies with halo 
masses less than 10 10 M have escape fractions enhanced by 
3-6 x over models without runaway stars. Various other mod- 
els find that the escape fraction varies (either hi gher or lower) 
with halo mass due to a v ariety of effects (e.g . JGnedin et all 
120081: IWise & Cenl 120091). inc luding holes in the ISM from 
supernovae (e.g. jYaiima e t al. 2009J). Thus, a relatively high 
escape fraction, and an escape fraction that varies with halo 
mass and redshift appears justified, though it remains to be 
verified observationally. 

In the following discussion, we will use the value of C/f esc = 
10 as our canonical model, which represents a scenario of C = 
3 and f t , sc . = 30%, consistent with the above results. Should 
the reader prefer different values, we show on Figure 3 critical 
curves for sustaining a fully reionized IGM for four values of 
C/f eJC : 3, 10, 30 and 100. In our discussion, we will discuss 
the impact of a different escape fraction on our conclusions. 

4.2. Can the Observed Galaxies Reionize the Universe ? 

In Figure 3, we show three previous results at z = 6, all 
derived solely using A CS optical imaging. The first, from 
iGiavalisco et al.l ( 12004b is based on a sample of 1 17 i'-dropout 
galaxies in the two GOODS fields. Correcting their sample 
for incompleteness down to Mj/y = -19.3, they found puy = 
1.15 ± 0.23 x 10 26 erg s" 1 Hz -1 Mpc" 3 . To convert this to 
our common magnitude limit of Mf/y = -18, we use the z = 
6 luminosity function of lBouwens et alj ( 120081) . finding a cor- 
rection factor of 2.0 (i.e., p uv [M < -18] = 2.3 ± 0.46 x 10 26 
erg s -1 Hz -1 Mpc -3 ). We note that in their study, the exist- 
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FIG. 3. — The specific luminosity density (puv) versus redshift, where our observed samples, corrected for incomplet eness down to Misoo < —18, are plotted 
as the cyan ci rcles. Results from the literature are shown as squares at z — 4, 5 and 6 from th e GOODS fields (red; [Giavalisco et al. 2004]), z = 6 from the 
HUD F (green IBunker et al.ll2004h z = 7 and 8 from the HUDF (purple; IFinkelstein eTaI]|2fjTo|) and z = 6 from GOODS+HUDF and z = 10 from the HUDF 
(blue; IBouwens et a l. 2006, 201 la]). All squares have been adjusted to a limiting magnitude of —18 (with the exception of the z = 7 and 8 point s, which were not 
corrected for incomple teness). The inverted triangles denote the in tegrated luminosity func tions down to M/,„, = -13 of[Bouwens et al. 12007]) at z = 4, 5 and 6, 
Bouwens et al. 1 20 1 lc) at z = 7, Bradley et al. 12012) at z = 8, and Bouwens et al. (201 la) at z = 10, where the gray error bar denotes the uncertainty on these 
integrated values due to the uncertainty in the Schechter function parameters (particularly the faint-end slope a in the higher redshift bins). The wide gray curves 
denote the v alue of puv needed to sustain a fully reionized IGM at a given redshift, for a given ratio of the clumping factor C over the escape fraction of ionizing 
photons f esc iMadau et al. 1999). The width of these curves represent changing stellar metallicities, from 0.02 < Z/Zq< 1.0. The right-hand vertical axis shows 
the corresponding intrinsic number density (i.e., prior to escape) of ionizing photons for a given specific UV luminosity density, using the median of a range 
of ages and metallicities, and assuming a constant star-formation history. This axis can be multiplied by the reader's choice o f f esc . The green curve shows the 
predicted luminosity density for sources brighter than Muoo = — 18 from the hydrodynamic simulations of Finlator et al. (201 1). 



ing depth in the z' data was 26.5; the data used in our present 
study is a full magnitude deepen 

The second observation, from lBunker et al. (2004) is based 
on a sample of 54 i'-dropout galaxies in the HUDF, finding 
p uv [M < -18.2] = 0.357 ± 0.071 x 10 26 erg s" 1 Hz" 1 Mpc" 3 ; 
we apply a factor of 1 .08 to convert this to our common mag- 
nitude limit. Finally, the third study, from iBouwens et al.l 
J2006T) uses both of the aforementioned datasets, including 
also the HUDF parallel fields, for a total sample of ~500 
galaxies, finding p uv [M < -17.5] = 1.77 ± 0.45 x 10 26 erg 
s" 1 Hz" 1 Mpc" 3 ; we apply a correction factor of 0.85 to these 
observations to convert to Pw\M < -181 . 

Interestingly, while the IBunker et al.l (120041) survey goes 
deeper, d own to Mpy ~ —18.2, t han the shallower GOODS 
data from Gi avalisco et al.l ( 120041) . one would have expected 
the former to yield a larger value of puv (prior to the cor- 
rection), but the opposite is observed. After we correct all 
three studies to a common magnitude limit, we would expect 
them to be in agreement, which is not the case, as they differ 
by more than a factor of five. There are a few reasons why 



such a discrepancy i s not surprising. The first is due to cos- 
mic variance, as the IBunker et al. (2004]) study only probed 
the 1 1 arc min 2 ACS data in the H UDF, while lGiava lisco et al. 
(120041) and IBouwens etalj (120061) studies included the ~ 300 
arcmin 2 GOODS fields. Secondly, as these galaxies were se- 
lected only with ACS imaging in all three surveys, they were 
observed in only a single band (z')- This can be dangerous, as 
the likelihood of selecting a noise peak is greatly enhanced 
when using a single b and detection versu s dual band (see 
IBouwens et alj d201 lal) and lYan et alj d2010h for a similar dis- 
crepancy in single-band selected z = 10 galaxy number counts; 
though the z = 10 disconnect is greatly enhanced due to the 
extreme faintness of such sources). 

The time is thus ripe to re-investigate the luminosity den- 
sity at z = 6, using our up-to-date dataset. Our current sam- 
ple of z = 6 galaxies is more robust, as our selection over 
both the GOODS-S field and the HUDF09 fields covers both 
a large volume, and a large dynamic range in object bright- 
ness. Additionally, the availability of the WFC3 data allows 
us to observe z = 6 candidate galaxies in four bands, increas- 
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ing the robustness of the sample. In particular, observations 
at four wavelengths throughout the rest-frame UV yield much 
more robust estimates of Myy than the single point previ- 
ously available at z = 6. We find 3 02 z = 6 galaxy candi- 
dat es, or nearly double t he size of the iGiavalisco et al.l (120041) 
andlBunker et al.l (120041) samples combined (though less than 
iBouwens et al.l (120061) : we expect this to be rectified when 
we soon incorporate the incoming CANDELS data in the 
GOODS-N field). We measure p uv = 0.79 ± 0.06 x 10 26 erg 
s" 1 Hz" 1 Mpc~ 3 , which is in between the previously measured 
values discussed above, with much smaller uncertainties. 

As shown in Figure 3, the observable (i.e., M1500 < -18) 
z = 6 galaxy population is capable of sustaining reioniza- 
tion for our canonical model of a low clumping factor and 
a moderate average escape fraction (C = 3 and f esc = 30%). 
Thus, in a change from a number of prev ious studies (e.g., 
lYan & Windhorstll2004l:lBunker et al.ll2004l) . it is not required 
to invoke galaxies fainter than our detection threshold to com- 
plete reionization by z = 6, due to both our more robust 
measurement, as well as the updated assumptions about the 
clumping factor. Nor do we need to invoke an escape frac- 
tion of unity, as f esc ~ 30% on average is sufficient (if C « 3). 
Finally, we reiterate this is also independent of assumptions 
about the form and faint-end cutoff of the luminosity func- 
tion. We note that our inclusion of deeper z' imaging, multiple 
fields, as well as near-infrared imaging yields a more robust 
estimate of the luminosity density, in contrast to the previous 
observations, which differed by a factor of > 5 x . 

Different conclusions are reached at z = 7-8. As shown in 
Figure 3, the specific UV luminosity density experiences a 
steady decline with redshift from z = 4 - 8. Combined with 
the denser IGM at higher redshift, necessitating a higher puy 
to compensate for the shorter recombination time, the observ- 
able galaxy population at z = 7-8 falls short of what is needed 
to sustain a fully reionized IGM, unless the escape fraction is 
unity (which is unlikely to be the case for all galaxies). Asim- 
ilar result is seen in the simulations of lSalvaterra et al.l (1201 11) . 
as shown in their Figure 8. They find that galaxies observable 
with HST can sustain a reionized IGM if f esc = 20% and C = 
5 at z = 6, with C < 2 (or a higher escape fraction) necessary 
at z > 7 to maintain reionization with the observable galaxies. 
However, their simulation also probes down to Mj/y = -15, 
and at that level they find that the observable galaxies with 
HST only account for ~20% of the total number of escaping 
ionizing photons. 

At z = 7, we find that the observed galaxies only account 
for ~50% of the required UV luminosity density inferred to 
sustain reionization. As we discuss later in §4.7, observations 
reaching ~ 1 magnitude fainter may yield enough galaxies 
(assuming the published luminosity functions) to account for 
a fully reionized IGM if C/f esc = 10, sim ilar to conclusions at 
this redshift from lRobertson et al.l (120101 ). We explore the pos- 
sible contribution from galaxies below our detection threshold 
in §4.4, and examine possible changing escape fractions in 
§4.5. 

Our observations from only the observed galaxy population 
are consistent with independent evidence that the universe 
was not comp letely reionized a t z > 7 . Studies of both Lya 
emission (e.g-lStark et al.l l2010Ll20Tll) and quasar near-zones 
(e.g., iFan et al. 1 120061) at~z ~ 6 yield results consistent with 
a negligible neutral fraction, implying that reionization was 
complete by z = 6. However, recent spectroscopic followup of 
z ~ 7 galaxies has yielded a much lower success rate of Lya 



emission detection than there should have bee n 
EW distributions at z = 4 - 6 dStark et al.ll20"T 



I N201 It lOno et al.H2012USchenker et al. 
2012 ; lTr*eu et al.H2012 ). implying neutral fractions as high as 



riven the Lya 
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~50%. A similar decrease in the Lya luminos ity function 
has b e en measured fr om z = 5.7 to z = 6.6 (e.g.jOuchi et al. 
I2010t IHu etail 120101) , though see also iMalhotra & Rhoads 
(120041) . Followup high-resolution spectroscopic study of the 
near-zone ar ound the recently disc overed z = 7 quasar ULAS 
Jl 120+0641 dMortlock et al.ll201 11) yields evidence for a neu- 
tral fraction > 10%. We note that the Lya results are also de- 
pendent on the neutral gas distribu tion within the galaxy (e.g., 
see §6.3 of Finkelstein et al. 20jTj), an d the quasar results are 
based on a single sightline , which could be bi ased by a nearby 
cloud of neutral gas (e.g., B olton et al.ll2011l) . thus more data 
are needed to conclusively show that the neutral fraction is in- 
deed rising. In any case, if our assumptions of a low clumping 
factor and a moderate escape fraction are reasonably correct, 
than the combination of galaxy, quasar and spectroscopic ob- 
servations point to an end to reionization by z = 6, but not 
earlier than z ~ 6.5. This then implies that through some 
combination of limiting magnitude and luminosity dependent 
escape fractions, the observed galaxy populations are an ex- 
cellent tracer of the ionizing photon background. 

In order to place these results in context, we e xamine the 
recent study of lKuhlen & Faucher-Giguerel ( 120121 ). who com- 
bined the WMAP and SPT results with the observed luminos- 
ity function evolution (extrapolated to higher redshift), and 
z < 6 Lya forest measurements to investigate a concordance 
model of reionization. They find that the available data pre- 
dict that the universe was ~50% ionized at z ~ 10. They note 
that, given the (less co nservative) Az < 4.4 constraint from 
the SPT measurements dZahn et al.ll2.01 ll) . scenarios invoking 
a very high contribution from very faint galaxies will result in 
a reioniz ation history which is too extended. However, recent 
work by dMesinger et al.| [2012) shows that there likely do ex- 
ist correlations between the thermal SZ power and the cosmic 
infrared background; in the presence of these correlations, the 
con straint on the duration of reionization degrades to Az < 
7.9 (IZahn et al.l201 UlMesinger et al.l20l2h . As shown in Fig- 
ure 3, if the IGM is truly ~ 50 ionized at z ~ 10, then there 
will need to be a significant contribution from faint galaxies, 
necessitating a longer duration of reionization. These conclu- 
sions thus stand to be revisited when a better understanding 
of this effect and its possible dependance on correlations are 
obtained. 

We note that our quoted re st-UV luminosities are not cor- 
rected for dust attenuation. In lFinkelstein et all (1201 ll) . we ex- 
amined the average dust attenuation as a function of redshift, 
finding that z = 6 galaxies appeared to contain little dust, and 
z = 7 galaxies appeared to be dust free (due to only two de- 
tected fluxes, robust estimates of the extinction at z = 8 were 
difficult, though they were consistent with little-to-no dust). 
Using the values of the UV spectral slope /3 from that work, 
and the conversion from /? to the attenuation at 1600 A from 
iMeurer etafl (119991 see discussion in Finkelstein et al. 201 1 
for the applicability of this relation at high redshift), we find 
that A1600 = 0.3 1 mag at z = 6, and A1600 = mag at z = 7 and 
8. Thus, accounting for dust does not change the above con- 
clusions at z = 7 and 8, while the observed specific luminosity 
density may in fact be up to 30% higher at z = 6. However, 
any dust present will likely also attenuate the ionizing pho- 
tons, though the dust attenuation curves are not well sampled 
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at such blue wavelengths lCalzetti et al.l (e.g. . 11994b . In light of 
this, we conclude that using the observed luminosity density 
is appropriate, assuming that the ionizing radiation is not sig- 
nificantly more (or less) attenuated than the rest-frame 1500 
A light. 

Finally, in Figure 3 we show in green the predicted luminos- 
ity density for galaxies at M1500 < -18, derived from a radia- 
tion hydrodynamic simulation that incorporates both galactic 
outfl ows and a self - consis tently grown ionizing background 
as in iFinlator et al.1 (1201 ll) . while subtending a 12 h" 1 Mpc 
volume at the same resolution. At all redshifts, these predicted 
values fall short of our observed luminosity densities. Some 
of this discrepancy may be due to volume limitations, but as 
we show below, the sub-L* galaxies are the dominant con- 
tributors, hence they should be reasonably well-represented. 
This discrepancy could indicate that a higher star-formation 
efficiency in the models is warranted, perhaps via weaker out- 
flows. It also relaxes the requir ements for the ionizi ng es- 
cape fraction. In particular, the IFinlator et al.l (1201 lb mod- 
els required an escape fraction of 50% in order to complete 
reionization by z = 6; weakening feedback by a factor of 2-3 
would suppress the required escape fraction by the same fac- 
tor, bringing it closer to the observed values at lower reds hift 
dShaplev -et al.ll200^[sian7et^lMllNestor et al.ll20"TTl) . as 
well as those we find below in §4.5. 

4.3. Contribution of Different Luminosity Ranges 

We now investigate how the different luminosity ranges of 
our observed galaxy sample contribute to the total specific UV 
luminosity density. We reiterate that this is independent of 
any luminosity function parameterization (with the exception 
of the definition of the characteristic luminosity L*). In Figure 
4, we show the fraction of the critical value of pyy necessary 
for sustaining a fully reionized IGM reached by our observed 
population at each redshift, for all galaxies (large circles), and 
for galaxies with luminosities greater than 0.2, 0.5 and 1 x L* 
(small colored symbols), where the appropriate value of L* 
for each redshift is used. 

At z = 6, this figure reinforces our observation that the ob- 
servable population at z = 6 appears sufficient to maintain an 
ionized IGM, while those at z = 7 and 8 are not. At z = 6, the 
ratio of the specific UV luminosity density from galaxies with 
L > 0.2L* to that of galaxies with L > L* is a factor of —2.6. 
This highlights that it is very much the sub-L* galaxies which 
are dominating the reionization photon budget. Additionally, 
galaxies with 0.2L* < L < 0.5L* contribute a roughly equal 
amount to reionization as those with L > L* . These ratios are 
similar at z = 7. 

At z = 8, the results are even more tilted in favor of the 
fainter galaxies, with the ratio of puv from galaxies with L 
> 0.2L* to that of galaxies with L > L* at a factor of —4.8, 
and galaxies at 0.2L* < L < 0.5L* contributing double the 
number of UV (and plausibly ionizing) photons as the L > L* 
galaxies (though both woefully short of the critical value of 
Puv)- Our observations thus imply that the faint-end slope of 
the luminosity function at z > 6 is steep, and that the slope 
is becoming even steeper as z — > 8 (in agreement with pre- 
vious studies; e.g., IBouwens et alj|201 lbtlOesch et al1l2012t 
iBradlev et al.ll2012T) . 

4.4. Galaxies Below our Detection Threshold 

Although our observed population appears consistent with 
multiple lines of evidence pointing to the completion of reion- 
ization by, but not too long before, z = 6, our conclusions 
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FIG. 4. — The fraction of the critical specific UV luminosity density nec- 
essary for sustaining reionization (denoted by the gray bar for C/f esc = 10) 
obtained from the (incompleteness-corrected) observed galaxy population 
(large cyan points) as a function of redshift, assuming C/f eS c = 10 (left ver- 
tical axis) and 30 (right vertical axis). The blue inverted triangles, green 
squares, and red triangles show the critical fraction for galaxies with lumi- 
nosities greater than 0.2, 0.5 and 1 X the characteristic luminosity L* , re- 
spectively. At z = 6, where the observed population can reionize the universe, 
this figure highlights that while galaxies below our detection limit are not nec- 
essary, it is the fainter portion of the observed population which is dominating 
puv, as galaxies with luminosities greater than 0.2L* contribute ~2.6x more 
than L > L* galaxies. This is even more so the case at z = 7 and 8, though at 
these redshifts the entire observed population is short by a factor of ~2.5-3 
of the ionizing photons necessary to sustain reionization for C/f esc = 10. 

are dominated by assumptions about the escape fraction of 
ionizing photons. Thus, it may be that fainter galaxies do 
play a large role. Although we will not observe them until 
the launch of the James Webb Space Telescope (JWST), it is 
highly likely that there do exist galaxies below our observa- 
tional limit. However, this introduces another uncertainty, in 
that we do not know how far down the luminosity function 
galaxies are able to form. 

At some limiting dark matter halo mass, the gas within a 
halo will no longer be able to self-shield against an ioniz- 
ing background from the sources which began reionization. 
In such halos star formation will be suppressed, effectively 
cutting off the UV luminosity function at a limiting magni- 
tude Mn,,, (e.g .,| Shapiro et al.l|1994tlThoul & Weinberg|[T996 : 
| Gnedirll200(l IPiikstra et al.l I2004bt lOkamoto etaU 12008 ; 
IFinlator et al.l 1201 11) . though in the case of inhomogeneous 
reionization, very small halos may still contribute a signif- 
icant number of ion izing photons before they are quenched 
( IFinlator et al.l 120 111) . This has not yet been seen observa- 
tionally, as there is no evidence for a turnover i n the lumi- 
nosity function (e.g.. IBouwens et aTll2007l f20 1 Id) , which im- 
plies that such a cutoff occurs at M/; m > -16.5 (the limiting 
absolute magnitude for z = 4 galaxies in the HUDF). Such a 
turnover has also not been seen in l arge cosmolog i cal sim ula- 
tions. For example, the models of IFinlator et al.l (1201 lb find 
an increasing luminosity function down to their resolution 
limit, which corresponds to M/,„, = -13. Interestingly, a recent 
st udy of high redshif t gamma ray burst (GRB) host galaxies 
bv lTrenti et al.l (12012b finds that the non-detection of five z > 5 
GRB hosts implies that galaxie s fainter than Muy = -15 were 
likely present at z > 5 (see also lTanvir et al. 2012|). 

The inverted triangles in Figure 3 show the specific 
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UV luminosity density if we take the luminosity functions 
from the literature, and integrate them down to M ;, m = -13 
dBouwens et al.ll2007ll20TTal5: iBradlev et al.ll2012l) . The un- 
certainties on these values were derived via 10 3 Monte Carlo 
simulations, where in each simulation each Schechter func- 
tion parameter was modified by a Gaussian random number 
multiplied by the associated parameter uncertainty. 

We note that in these simulations, we have used the pub- 
lished uncertainties on each Schechter function parameter; 
thus, they were free to vary independent of each other. How- 
ever, these parameters, particularly the characteristic magni- 
tude M* , and the faint-end slope a, are correlated, thus the 
uncertainty on the integrated luminosity functions may be 
smaller than those derived here. If the published luminosity 
functions accurately represent the true galaxy population, and 
the population does indeed extend this faint, than the existing 
galaxy population at z = 7 and z = 8 can sustain reionization 
for C/f esc =10. However, both luminosity functions are highly 
uncertain at this point, especially at z = 8, where the error bar 
stretches well below the value necessary to sustain reioniza- 
tion. These luminosity function uncertainties, combined with 
the uncertainty on the value of M/„„, render the contribution 
of the unseen faint galaxies to reionization difficult to assess. 

4.5. Constraints on the Escape Fraction 

As discussed above, if the escape fraction is modestly high, 
at ~30%, then the observable galaxy population at z = 6 can 
sustain a fully reionized IGM. However, it is highly probable 
that galaxies significantly fainter than this limit exist, given 
results from a wide variety of cosmological simulations (e.g. , 
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20T0t . ISalvaterra et al.ll201ll IFinlator et al.H201 It 
120 12|). and plausible , though scant, observational 



evidence (e.g. jTrenti et al. 20T2h. 

In light of the likelihood that the galaxy luminosity function 
extends fainter than we can observe, it is instructive to exam- 
ine the effect of our choice of 30% for the escape fraction in 
the above discussion. As shown in Figure 3, our observed 
galaxy population will yield just enough photons to complete 
reionization at z = 6 if the escape fraction is 30%. If the lumi- 
nosity function extends significantly fainter than Mf/y = -18, 
and the escape fraction is on average 30% for all galaxies, 
then more ionizing photons will be produced than are neces- 
sary to complete reionization. In principle, this is not a prob- 
lem, as once the IGM is reionized, it will stay ionized as long 
as the critical number of ionizing photons are being produced. 

However, there is one observational constraint we have not 
yet strongly considered, which is the measurement of the ion- 
izi ng emissivity from observ ations of the Lya forest opac- 
ity. iBolton & Haehnehl d2007l) studied this at z = 6 combining 
simulations with high-resolution spectroscopic observations 
of z > 6 quasars, and found that only 1.5-3 ionizing photons 
per hydrogen atom were created, leading them to believe that 
reionization was "photon-starved." At z = 6, they calculated 
a la upper limit on the emission rate of ionizing photons per 
comoving Mpc of < 2.6 x 10 50 s" 1 Mpc~ 3 . This result, though 
based on a relatively small number of sight-lines, implies an 
upper limit to the number of ionizing photons created at z = 6. 

In Figure 5, we compare this constraint from the Lya for- 
est to our observations, for a variety of escape fractions. We 
convert our observed UV luminosity density into a comoving 
rate of ionizing photons using the same set of stellar popula- 
tion model assumptions used to derive 653 in §3.4. Assum- 
ing Z = 0.2Z©, we find N io „ = 1.9 x 10 25 x pw x fesc- This 



assumes a star-formation history that is in equilibrium (i.e., 
massive O stars die as fast as they are formed). We use this to 
show in Figure 5 the inferred emission rate of ionizing pho- 
tons from both our ob served galaxies (in cya n), and the z = 6 
luminosity function of lBouwens et al.l (12008) integrated down 
to M[/v = -10 and -15 (light gray band, where the left and 
right edges of the band denote M/,- m = -10 and -15, respec- 
tively). 

If M/,„, ~ -18, an escape fraction of 30% is allowed within 
the 2a upper limit on the Lya forest constraint. With this 
constraint, the 2a (la) upper limit on the escape fraction is < 
34% (17%). As shown by the light gray curve, if the lumi- 
nosity function extends fainter than our observations, then the 
escape fraction must decrease, otherwise the Lya forest con- 
straints will be violated. Using the limit of Mw = -13 used 
above in §4.4, we find that accounting for galaxies fainter than 
our detection limit the average escape fraction must be < 13% 
(6%) at 2a (la). We note that the 2a limit on f esc in both cases 
(observed and integrated luminosity function) are consistent 
with sustaining reionization if the clumping factor C = 3. If 
future observations constrain the ionizing photon production 
rate to be closer to the la limit, then the galaxy population 
can only sustain reionization if C < 2. 

Given the current observations, we can thus say with 95% 
confidence that the average escape fraction is < 34%. Fu- 
ture observations with JWST will push the luminosity func- 
tions several magnitudes lower, and thus put more robust con- 
straints on this quantity. We however note that current obser- 
vations constraining the ionizing background at z = 6 are due 
to a relative low number of individual sight-lines to distant 
quasars. Additionally, these quasars reside in biased regions, 
thus their spectra may not represent a true probe of the mean 
IGM neutral fraction. Ongoing and future studies such as 
UKIDDS and VIKING will yield much greater probes of the 
ionizing background at z = 6-7, which can then be turned into 
stronger constraints on the escape fraction of ionizing photons 
from the z = 6 galaxy population. 

Finally, while the above discussion assumes the escape 
fraction represents an average of all galaxies, this does not 
have to be the case. It is more accurately a luminosity- 
weighted escape fraction; essentially, this quantity tells us 
what fraction of the total number of ionizing photons pro- 
duced in all galaxies escape into the IGM. It may also be the 
case that the escape fraction is time dependent; ionizing pho- 
tons may only escape from galaxies due to specific temporal 
events, such as interactions or starbursts. Our constraints limit 
the total number of ionizing photons which escape from all 
galaxies, which can be consistent with a time, luminosity or 
mass dependance, or no dependance at all. 



4.6. Constraints on the Reionization History 

As we move into the reionization epoch, and the ionized 
fraction in the universe becomes less than unity, one can use 
Equation 2 to calculate the luminosity density necessary to 
maintain a given ionized fraction by varying the volume ion- 
ized fraction xhh- We can thus indirectly estimate xhh as a 
function of redshift by dividing our observations by the criti- 
cal value of puv necessary to sustain a fully ionized IGM as a 
function of redshift. Thi s approach is a valid approximation to 
semi-analytic modeling dMadau et al Jl 19991) because the IGM 
recombination time remains less than the Hubble time for C = 
3 at z > 6. This is subject to the same assumptions as the 
above analysis, including: 
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FIG. 5. — Left) The inferred emission rate of ionizing photons from our (incomple teness-corrected) obse rved sample of z = 6 galaxies (cyan; the width of curve 
denotes the 68% confidence range), as well as from the z = 6 luminosity function ofBouwens et al. 1 2003) integrated down to Muv = ~ 10 (left side of light-gray 
curve) and -15 (right side). Both quantities are plotted as a function of the escape fraction of ionizing photons. The red region deno t es the constraints on this 
quantity at z = 6 from the Lya forest, with the lower and upper edges representing the 1 and 2a constraints from Bolton & Haehnelt (2007), respectively. The 
dark gray bands denote the emission rate of ionizing photons necessary to sustain reionization for a given value of the clumping factor. Right) The striped curve 
shows the allowed average escape fraction of ionizing photons from galaxies as a function of the limiting UV absolute magnitude, assuming the z = 6 luminosity 
function of Bouwens et al. 1 2007'), based on the 2a constraints from the Lya forest. The cyan star shows the allowed escape fraction from our observations with 
its associated uncertainty. If the galaxy luminosity function turns over shortly faintward of our observations, then the average escape fraction of ionizing photons 
is constrained to be f esc < 34% (17%) at 2<r (la). If the luminosity function extends down to Muv = _ 13, then this constraint tightens to f esc < 13% (6%) at 2a 
(1<t). These escape fractions can still reionize the IGM if the clumping factor is ~ 2-3. 



• Galaxies are the only source of ionizing photons in 
the z > 6 universe. While other sources (i.e., the rare 
quasar) may exist, their overall contribution is likely 
negligible. 

• The faint-end cutoff of the luminosity function. As 
above, we will consider two scenarios: 1) Only galax- 
ies above the observational limit of Muv = - 18, and 2) 
The published luminosity functions integrated down to 
M uv =-13. 

• The clumping factor of ionized hydrogen, which we 
take to be C = 3. 

• The escape fraction of ionizing photons; per §4.5, we 
assume here f esc = 30% when considering only the ob- 
servable galaxies, and f esc = 10% when considering 
the integrated luminosity functions. Both scenarios are 
consistent with the Lya forest constraints, and both can 
sustain a fully reionized IGM at z = 6 if C < 3. 

• We assume that the ratio C/f esc does not evolve over the 
redshifts we study. 

We plot the results of this analysis in Figure 6. Our 
incompleteness-corrected measurement from the observed 
galaxies (assuming f esc = 30%) is consistent with a volume 
ionized fraction of 1.0 at z = 6, ^0.4 at z = 7, and ~0.2 at 
z = 8. When we consider the integrated luminosity function, 
with f esc = 10%, we find a much wider range of possible vol- 
ume ionized fractions. We note that while the UV luminos- 
ity density from the integrated luminosity function will al- 
ways be larger than that for galaxies above our observational 
limit, here we assume two different escape fractions for the 



two magnitude limits, thus their ionizing luminosity densities 
are comparable. At z = 6, the integrated luminosity functions 
are consistent with ionized fractions of ~ 0.7 - 1.0. At z = 
7 and 8, this is reduced to 0.3 - 1.0, and 0.1 - 0.95, respec- 
tively. The very wide range of possible ionized fractions from 
the integrated luminosity functions are due to the increasingly 
larger uncertainties in the Schechter function parameters as 
one pushes to higher redshift. 

Also on Figure 6, we plot a number of results from the lit- 
erature, including con straints on the v olume ion ized fraction 
from quasars at z < 6 dFan et alj|2006l) and z = 7 dBolton et alj 
1201 II) . W MAP dKomatsu et al JI20 1 11) . Lya emitter luminosity 
functions dOuchi et al.ll2010[). and Lya spectroscopic studie s 
dPentericci et al.l201 UlSchenker et al.l2012l;IOno et al.l2012l) . 
Both our observed galaxy population, and the integrated lu- 
minosity functions, are consistent with this wide and com- 
plementary range of reionization probes. The only exception 
is the single known quasar at z = 7, which has a near-zone 
size consistent with an ionized fraction of < 90%, and ap- 
pears slightly inconsistent with our observed galaxy popula- 
tion (though consistent with the integrated luminosity func- 
tions). However, this is but a single sightline to the only cur- 
rently known z > 7 quasar, thus many more sight-lines are 
needed to truly constrain the IGM from this type of probe. 

Our new observations combined with the previous, comple- 
mentary, probes of the ionization state of the IGM show a con- 
sistent reionization picture, where the ionized fraction in the 
IGM is essentially unity at z < 6. At z ~ 7, many lines of evi- 
dence indicate a non-negligible neutral fraction, including our 
observed galaxy population. The integrated luminosity func- 
tion prefers an ionized fraction < unity, but it is consistent 
with a fully ionized IGM at la. This is not the case at z = 8, 
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FIG. 6. — The volume ionized fraction, xhii, of the IGM which can be sustained given the observed UV luminosity densities (and our assumptions on C and 
f eS c) as a function of redshift. The cyan and gray curves denote the value of xhii inferred from comparing our observed galaxies, and the integrated luminosity 
functions, respectively, to the critical luminosity density from Equation 2. We have assumed C/f esc = 10 when only considering the observable galaxies, and 
C/fesc = 30 when considering the integrated luminosity functions (both of which are consistent with limits from the z - 6 Lyct forest presented in §4.5, for C = 
3). We note that while the luminosity density from the integrated luminosity functions will always be larger than that for galaxies above our observational limit, 
here the two different escape fractions for the two magnitude limits result in their ionizing luminosity densities being comparable. We assume that galax ies are 
the only si gnificant sources of ionizing photons in this epoch. We plot constraints on xhii from spectroscopy of quasars at z < 6 from Fan et al. (2006) and at 
z = 7 from Bolton et al. 1 201 1) in red. The blue circle denotes constraints on xhii from the evolution in the Lya luminosity function from z = 5.7 to 6 . 6 from 
lOucMet al. (2010), w hile the blue bar d enotes the range of xhii values inferred from the z = 7 follow-up Lya spectroscopic studies of Pentericci et al. (201 1); 
ISchenker et al . (2012); Ono et al. |2012|)- The instantaneous redshift for reioniza tion from WMAP (10.6 ± 1.2) is indicated by the orange rectangle. The derived 
50% and 90% xhii redshifts from the study of Kuhlen & Faucher-Giguere (2012) are shown in green, where the width of the green bars denote the lcr uncertainty 
(from Figure 9 of Kuhlen et al.). Our observations, which are consistent with a wide range of other, complementary, approches, find a picture where the universe 
is fully ionized by z = 6, with the neutral fraction becoming non-negligible at z > 7. This is independent of whether we consider only the observable galaxies, or 
the integrated luminosity functions, each with their own permitted escape fraction. These observations are in moderate tension with the concordance models of 
Kuhlen et al. at z > 9, as these models prefer an earlier redshift for xhii = 50%. This may be resolved with future deep observations, which will provide much 
better constraints on the z > 8 luminosity functions, and show if the strong drop-off in luminosity density at z > 8, currently inferred from a single source, is 
accurate. 



where galaxies are the only probe we currently have, and they 
seem to indicate an ionized fraction of < 80-90% (unless the 
escape fraction in all galaxies is unity, which would indicate 
extreme evolution from z = 6, and is unlikely; see §4.5). Un- 
fortunately, more robust constraints on the evolution of X//// 
from the integrated luminosity functions are difficult due to 
the uncertainty in the Schechter function parameters. We note 
that we have used the uncertainties on each Schechter func- 
tion parameter to compute the uncertainty in the integrated 
luminosity function (see §4.4 for more details). These param- 
eters, particularly the characteristic magnitude M*, and the 
faint-end slope a, are known to be correlated, thus the uncer- 
tainty on the integrated luminosity functions may be smaller 
than those derived here. In any case, at z > 7, the faint-end 
slope is not well constrained, which is responsible for the bulk 
of the uncertainty on X////. More robust luminosity functions 
in this epoch will help decrease the uncertainty on Xhh from 



galaxy measurements. 

Finally on Figure 6, we also show the concordance models 
o flKuhlen & Faucher-Giguerel(l2012l) . which present a picture 
of the reionization history of the universe consistent with the 
Lya forest, WMAP optical depth, and the kSZ effect. They 
also fold in the observed galaxy luminosity functions, though 
in this case they assume that the luminosity function parame- 
ters evolve smoothly with redshift; we do not know if this is 
truly the case, particularly at z > 8, beyond which the strong 
drop-off in the observed luminosity density is intriguing. We 
plot their results for the redshift of a 50% and 90% ionized 
fraction. Comparing our inferred values of X//// from galaxies 
with their conclusions, we find general agreement at z < 8, in 
that they find that reionization completes by z = 6, and that it 
is 90% complete at z = 7.5 ± 1.5. 

However, at higher redshift, there appears to be some ten- 
sion between our observations and their results, in that the 



14 



Finkelstein et al. 




C/fesc 

FIG. 7. — The colored curves show the limiting UV absolute magnitude for different redshifts necessary to sustain reionization for a given ratio between the 
IGM clumping factor and the ionizing escape fraction. The hatched region denotes the magnitude range where star formation is likely to be suppressed via 
photo-heating from the sources which began the process of reionization at higher redshift. The dashed line shows our canonical model of a clumping factor of 
3 and f esc = 30%, or C/f« c = 10, which is a reasonable assumption for the conditions in the universe at z > 6. The horizontal gray bar denotes the limiting UV 
magnitude of HST observations of the HUDF at z = 7 (m^ = 29). A hypothetical JWST deep field will reach 2.5 mag deeper, to m^g = 31.5, or M(/v = -15.5 (for 
Z = 7). At z = 6, our observations (shown in Figure 3) already reveal a population of galaxies capable of sustaining reionization. At z = 7, the observed population 
is not bright enough to sustain reionization, thus if reionization completes at z > 6, fainter galaxies are necessary, which should be revealed by JWST. At z = 8, 
even JWST cannot image deep enough to see the hypothetical extremely faint population necessary to complete reionization at z = 8. 



iKuhlen & Faucher^G iguere (2012) study prefers a 50% reion- 
ized universe at a higher redshift than appears possible with 
the observations, though it is marginally consistent with the 
integrated l uminosity functions. However, we stress that the 
results from lKuhlen & Faucher-Giguerel(l2012l) were tuned to 
match the observed luminosity functions. Any tension likely 
arises from the WMAP constraint on the Thomson scattering 
optical depth due to electrons, which we have not accounted 
for. However, as the WMAP constraint prefers a redshift of 
instantaneous reionization of 10.6 ± 1.2, it is not surprising 
that observations showing a strong decline in the galaxy UV 
luminosity density at z > 8 would underproduce the electron 
optical depth. This can be reconciled by allowing the escape 
fraction, and/or the limiting magnitude to evo l ve wi th red- 
shift, as is done by IKuhlen & Faucher-Giguerd ( 120121) . This 
would imply that very faint galaxies may be necessary to drive 
reionization at the earliest stages, and that reionization gets 
its start at much higher redshifts (i.e., z > 10). We point out 
that the z = 8 luminosity function, and the limits on the z = 
10 luminosity function, are very rough at this point in time. 
Future observations, likely with /W5T(see §4.7), will greatly 
improve constraints on the luminosity functions in this early 
epoch, and at that time we can learn more about whether the 
contribution of galaxies to the early stages of reionization are 
consistent with all other constraints. 

4.7. Implications for Future Surveys 

At z = 6, our observations (shown in Figures 3 and 4) al- 
ready reveal a population of galaxies capable of sustaining 



reionization, while at z = 7 and 8 the observed galaxy popu- 
lation falls short. Our current observations are however sus- 
ceptible to the uncertainty about what lies below our detection 
threshold, as well as our assumptions about the escape frac- 
tion. The former can potentially be alleviated by integrating 
down the luminosity function, as we have done in §4.4, but 
this yields its own uncertainty: how faint does the UV lumi- 
nosity function continue? 

In order to measure this, we will require the deeper obser- 
vations only possible with JWST. In a hypothetical ~800 hr 
program, JWST can reach mAB(5cr) = 31.5 mag in four bands; 
F070W, F090W, Fl 15W and F150W; which can be used to se- 
lect samples at similar redshifts as in our study 3 . Is this deep 
enough to probe faint enough values of M/,„, to conclusively 
find the reionizing population at z = 7 and 8? To investigate 
this, we created Figure 7, where we plot the limiting absolute 
UV magnitude necessary to sustain reionization for a given 
ratio of C/f, „ r , assuming the lumin osity functions from the lit- 
erature (see Iwilkins et al J (120 1 1 ah for a similar analysis at z = 
7). 

In this figure, we show the limiting magnitude of our HST 
study, as well as that of the hypothetical JWST deep field, 
where each colored line denotes the value of M/„„ and C/f e „- 
necessary for a given redshift. As we discussed above, for 
our assumption of C/f esc =10, our z = 6 observed galaxy pop- 
ulation is sufficient to sustain a reionized IGM. In fact, this 
figure shows that the observed population can reionize the 

3 http://jwstetc.stsci.edu/etc/input/nircam/imaging/ 
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IGM for values of C/f esc < 14. At z = 6, JWST will trace 
the rest-frame UV luminosity function down to faint levels, 
enabling the search for a turnover at faint magnitudes down to 
M1500 = -15.3. If the luminosity function continues its steep 
climb to this magnitude, then as shown in Figure 5, the aver- 
age escape fraction for galaxies will be constrained to be less 
than -15%. 

At z = 7, galaxies brighter than HSTs magnitude limit can 
only reionize the universe if C/f esc < 6, which, for C > 3, 
necessitates an escape fraction > 50%. While this is pos- 
sible, there is not yet strong observational evidence to sup- 
port such a scenario. However, JWST will reach deep enough, 
down to M/,„, = -15.5; if the luminosity function continues 
its steep slope down to such faint magnitudes, and C/f esc < 
15, then galaxies could have r eionized the IGM by z = 7 
(see also (Salvaterra et al.l 1201 11) . At z = 8, even JWST can 
not probe deep enough to see the limiting magnitude neces- 
sary to reionize the IGM (again assuming C/f esc > 10). If 
the universe was reionized by z = 8, it necessitates a luminos- 
ity function that continues with a steep slope (a — -2) down 
to M/,„, = -12. However, such a result would create tension 
with the observed dearth of Lya emission, and the small size 
of quasar near zones at z = 7 (see §4.2). This tension could 
be resolved if the escape fraction of galaxies is luminosity 
dependent, and galaxies with L <C L* have very low escape 
fractions, w hich could be the ca se if they are highly gas dom- 
inated (e.g.. lGnedinetal.ll2008l) . 

5. CONCLUSIONS 

We have measured the specific UV luminosity density 
(Puv) from a large sample of galaxies at z = 6, 7 and 8. Using a 
theoretical model to convert from UV to ionizing photons, and 
determining the level of ionizing photons necessary to sustain 
a fully reionized IGM, we have investigated the contribution 
of observed galaxies to the reionization of the IGM. A large 
body of work over the past — decade points to a probable 
low clumping factor. Thus, previous surveys which assumed 
a large clumping factor and deduced that galaxies could not 
reionize the universe may have reached the opposite conclu- 
sion today. Unfortunately, not much is known about the es- 
cape fraction of ionizing photons at high redshift, but current 
observations are consistent with a rising escape fraction with 
increasing redshift, or at least an increased incidence of large 
escape fractions. With our new sample, and our increased un- 
derstanding of the process of reionization, we have reached 
the following conclusions: 

• We measure p uv at z = 6 to be 0.79 ± 0.06 x 10 26 erg 
s" 1 Hz" 1 Mpc" 3 for the observed population of galax- 
ies; this is much more robust than previous measures, 
which had much smaller samples and were limited to 
single band detections, and differed by a factor of > 5. 
Additionally, this result is independent of a Schechter 
function parametrization for the form of the luminos- 
ity function. Assuming C/f esc — 10, we find that the 
observed galaxy population at z = 6 is sufficient to sus- 
tain a fully reionized IGM. In contrast to many previ- 
ous studies, this result does not need to invoke galaxies 
fainter than our detection threshold. 

• We measure puy of the observed galaxy population at 
z = 7 and 8 to be 0.44 ± 0.04 and 0.27 ± 0.06 x 10 26 erg 
s" 1 Hz" 1 Mpc" 3 , respectively. These values fall short 
of what is needed to reionize the IGM. While fainter 



galaxies may render the population capable of complet- 
ing reionization, even JWST will not go deep enough to 
see the necessary limiting magnitude at z = 8. 

• If we only consider the observed galaxies, we find a 
scenario where though z = 7 and 8 galaxies contribute 
to reionization, the IGM only becomes fully reionized 
by z = 6. This is consistent with independent studies of 
both Lya emission from normal galaxies, and the study 
of Lya absorption in the near-zone of the z = 7 quasar, 
both of which find that the IGM at z = 7 likely has a 
neutral fraction > 10%. 

• Accounting for constraints on the total emission rate of 
ionizing photons from measurements of the Lya forest 
at z = 6, we find that if our observed galaxies repre- 
sent the total galaxy population, then the average z = 6 
escape fraction must be f esc < 34%. However, fainter 
galaxies likely do exist, and if the luminosity function 
extends down to M/,-„, = -13, then the average escape 
fraction constraint tightens to f esc < 13%. These es- 
cape fractions can sustain reionization z = 6 either with 
our observed galaxies, or with the integrated luminosity 
function, for values of the clumping factor < 3. 

• We investigate the evolution of the volume ionized frac- 
tion Xh// of the IGM. We do this by combining the ob- 
served galaxy population of our study and the lumi- 
nosity functions from the literature with limits on the 
escape fraction considering both limiting magnitudes. 
While the IGM appears to be fully ionized by z < 6, the 
volume ionized fraction appears to drop below unity by 
z = 7, consistent with a number of complementary anal- 
yses. Considering only the (roughly known) contribu- 
tion from galaxies at z > 8, the volume ionized fraction 
may drop substantially below unity, though this is in 
mild tension with reionization models which incorpo- 
rate the optical depth due to electron scattering. Bet- 
ter constraints on the z > 8 luminosity functions, likely 
from JWST, will place stronger constraints on the reion- 
ization history at that early time. If the strong drop in 
the luminosity density at z > 8 is verified, it may im- 
ply some additional source of ionizing photons at z > 
10 is warranted to remain consistent with the WMAP- 
measured Thomson scattering optical depth. 

We conclude that in order to make progress on the issue 
of when and how reionization happened, we need better con- 
straints on the escape fraction of ionizing photons or the limit- 
ing magnitude of the UV luminosity function. As these quan- 
tities are correlated due to the constraints from the Lya forest, 
constraints on one will yield more robust constraints on the 
other. While direct measurement of the escape fraction at high 
redshift is not possible, much deeper observations with JWST 
should confirm whether the luminosity function extends down 
to at least M/,-„, = -15.5, which will place stronger constraints 
on the escape fraction at z = 6, and uncover whether the z = 7 
galaxy population is sufficient to reionize the universe. 
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